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Abstract 

We investigate the scale-dependence of Eulerian volume averages of 
scalar functions on Riemannian three-manifolds. We propose a comple- 
mentary view of a Lagrangian scaling of variables as opposed to their 
Eulerian averaging on spatial domains. This program explains rigorously 
the origin of the Ricci deformation flow for the metric, a flow which, on 
heuristic grounds, has been already suggested as a possible candidate for 
averaging the initial data set for cosmological spacetimes. 
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To Ruggiero, always open to learning new things 

1 Introduction 

Averaged inhomogeneous cosmologies are in the forefront of interest, since cos- 
mological parameters like the rate of expansion or the mass density are to be 
considered as volume-averaged quantities and only these can be compared with 
observations 0,[§,[||,@,[§,[§,[0,|!]. For this reason the relevant param- 
eters arc intrinsically scale-dependent and one wishes to control this dependence 
without restricting the cosmological model by unphysical assumptions. Con- 
sider a three-dimensional manifold equipped with a Riemannian metric (£,<?) 
as a part of the initial data set for a cosmological spacetime. On such a hyper- 
surface we may select a simply-connected spatial region and evaluate certain 
average properties of the physical (scalar) variables on that domain such as, e.g., 
the volume-averaged mass density field, or the volume-averaged scalar curva- 
ture. These (covariant) average values are functionally dependent on position 
and shape of the chosen domain of averaging. Let us now be more specific and 
relate the averages to scaling properties of the physical variables. It is then 
natural to identify the domain with a geodesic ball centred on a given position, 
and introduce the spatial scale via the geodesic radius of the ball. We may 
consider the variation of this radius and so explore the scaling properties of 
average characteristics on the whole manifold. Doing so for every position we 
arrive at smoothed fields that capture the effective scale-dependent dynamics 
and so naturally provide the theoretical values of observables that are drawn 
from a limited survey volume of a cosmological spacetime. 

We may call this point of view Eulerian - and everybody would first think of 
this point of view - in the sense that the spatial manifold is explored passively 
by blowing up the geodesic balls covering larger and larger volumes. On the 
contrary, the key idea of our approach consists in demonstrating that a corre- 
sponding active averaging procedure can be devised which is Lagrangian: wc 
hold the ball at a fixed (Lagrangian) radius, and deform the dynamical variables 
inside the ball such that the deformation corresponds to the smoothing of the 
fields. The variation with scale of, e.g., the density field or the metric is then 
mirrored by one-parameter families of successively deformed density fields and 
metrics. We shall show rigorously that this deformation corresponds to a first 
variation of the metric along the Ricci tensor, known as the Ricci flow. Since 
this flow has received great attention in the mathematical literature - major 
contributions are due to Hamilton Q , - we shall so translate the averaging 
procedure into a well-studied deformation flow, much in the spirit of a renormal- 
ization group flow ||, 0, Q, 0|. This analysis, putting on rigorous grounds 
the previous, largely heuristic use of the Ricci flow in smoothing out relativis- 
tic spacetimes, is part of a larger program aimed to a full-fledged analysis of 
averaging and scaling in relativistic cosmology. Such a program features a mass- 
preserving Ricci flow and its variants in order to provide a subtle technique for 
modelling an effective dynamics of inhomogencities in the Universe. 
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2 Averaging and Scaling put into a Geometrical 
Perspective 

In order to characterize the correct conceptual framework for addressing av- 
eraging and scaling properties in relativistic cosmology, let us remark that in 
general relativity we basically have one scaling variable related to the unit of 
distance: we can express the unit of time in terms of the unit of distance using 
the speed of light. Similarly, we can express the unit of mass through the unit of 
distance by using Newton's gravitational constant. This remark implies that the 
scaling properties of Einstein's equations typically generate a mapping between 
distinct initial data sets, (in this sense we are basically dealing with a renor- 
malization group transformation). As we shall see below in detail, a rigorous 
characterization of an averaging procedure in relativistic cosmology is indeed 
strictly connected to the scaling geometry of the initial data set for Einstein's 
equations. 

In order to discuss the scaling properties and the averaging procedure asso- 
ciated with an admissible set of initial data for a cosmological spacetime (we 
think of a setting within the Arnowitt-Deser-Misner formulation of general rel- 
ativity), a basic issue is to characterize explicitly a scale-dependent averaging 
for the empirical mass distribution g. This will be the main theme of this paper. 

2.1 Matter seen at many scales 

In order to characterize the scale over which we are smoothing the empirical 
mass distribution g, we need to study the distribution g by looking at its average 
behavior on (geodesic) balls on (£,<?) with different centers and radii. The idea 
is to move from the function g : S — > M + to an associated function, defined 
on S x K+ , and which captures some aspect of the behavior of the given g on 
average, at different scales and locations. The simplest function of this type is 
provided by the local volume average 



where p € £ is a generic point, d[i g is the Riemannian volume element associated 
with (£,g), and B(p;r) denotes the geodesic ball at center p and radius r in 
(£,c/), i.e., 



where d g (p,q) denotes the distance, in (£,</), between the point p and q. Note 
that if 




(1) 



B{p;r) = {qe(J:,g) : d g (p, q) < r } 



(2) 



diam = sup {d g {p 1 q) : p,q 



(3) 
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denotes the diameter of (E,g), then asr-> diam, we get 

at any point peE. Conversely, if g : E — > M + is locally summable then 



lim (e>B(p;r) = e(p), (5) 

for almost all points p e E. The passage from g to (g) B t p . r ^ corresponds to 
replacing the position-dependent empirical distribution of matter in B(j>;r) by 
a locally uniform distribution (g) B ( p . r ) which is constant over the typical scale 
r. Note that for the total (material) mass contained in B(p; r), we get 

M(B(p; r)) = ( gdn g = Vol (B(p; r)) (g) B[p . r) , (6) 

JBfp-r) 



and 



lim M(B(p;r))=M(E)= f gd^ g , (7) 

where M(E) is the total (material) mass contained in (£,<?). Our expectations 
in (£>) B ( p . r ) are motivated by the fact that on Euclidean 3-space (K 3 ,<5 a fc), if Q 
is a bounded function, then the local average (g) S ( p . r ) is a Lipschitz function 
on R 3 x M+, i.e., 

(e> B (p;r) - <e>B( g;r ) <C d h (p,q), (8) 

where Co is a constant and dh(p, q) is the (hyperbolic!) distance between p and 
<j. Thus, the local averages (g) B / p . r \ do not oscillate too wildly as p and r vary, 
and the replacement of g by {(£ , ) B ( p . r )} P eR 3 indeed provides an averaging of the 
original matter distribution over the length scale r. It is not obvious that such 
a nice behavior carries over to the Riemannian manifold (E,g). The point is 
that, even if the local averages {(g) B ( p . r )} P GS provide a controllable device of 
smoothing the matter distribution at the given scale r, they still depend in a 
sensible way on the geometry of the typical ball B(p\ r) as we vary the averaging 
radius. In this connection we need to understand how, as we rescale the domain 
B(p; r), the local average (g) B r p . r \ depends on the underlying geometry of (S, g). 
The reasoning here is slightly delicate, so we go into a few details that require 
some geometric preliminaries. 

Let us denote by 

ex Pp : T p E — > E (9) 
(v, r) i — ► exp p (rv) 
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the exponential mapping at p G (£,3), i-e., the map which to the vector rv G 
T p E ~ R 3 associates the point exp p (rv) G E reached at "time" r G K + by the 
unique geodesic issued at p G E with unit speed u G § 2 (1) • Let r > be such 
that expp is defined on the Euclidean ball 

B E (0,r) = {yeT p Z~R 3 : \y\ < r} (10) 

and exp p : Be{Q,t) — ► B(p,r) C S is a diffcomorphism onto its image. The 
largest radius r for which this happens, as p varies in E, is the injectivity 
radius inju of (E, g). Recall that a curve c Ptq : M D / — > (E,<?) connecting 
the base point p with the point g is called a segment, if its length L(cp >q ) is 
such that L(cp^q) — d g (p,q), and if it is parametrized proportional to arclength. 
The domain in T p Y> over which geodesies issued from p are segments (i.e., are 
distance-realizing) is called the segment domain and is defined according to 



Seg(p) = TpH; exp p (ru) : [0, 1] — > S is a segment] . (11) 
Note that Seg{p) is a closed star-shaped subset of T p £, and E = exp p (Seg(p) 



The exponential map acts diffcomorphically over the interior Seg(p) of Seg(p), 
and Seg(p) — Seg(p), the cut locus of p in T p E, is a set whose image under exp p 
has zero measure in E. In particular, it follows that 

M(E)= / g(y)d^= / e(ex Pp )ex P ;(^ 3 ), (12) 

JT. JSeg{p) 

where exp*(du g ) denotes the pull-back of the Riemannian measure du g over 
Seg{p) C TpTi under the exponential mapping. 

Let B(p; ro) denote a given geodesic ball of radius tq < injisi, and for vector 
fields X, Y, and Z, in (B(p;r),g) let R{X,Y)Z = R a bcd X c Y d Z b be the corre- 
sponding curvature operator. Since r : (B(p;r),g) — > M is a distance function 
(i.e., I Vr| = 1), the geometry of B(p;r) can be described by the following set 
of equations 12 : 

(Vg r S)(X) + S 2 (X) = -R(X,d r )d r , (13) 
(L dr g)(X,Y) = 2g(S(X) ) Y), (14) 
V dr S = L dr S, (15) 

where d r — Vr is the gradient of r, La^g is the Lie derivative of the metric 
in the radial direction d r , and the shape tensor S = V 2 r is the Hessian of r. 
Such equations follow from the Gauss-Weingarten relations applied to study 
the r-constant slices U r = {exp p (rv) G E : r = const.} which are the images 
in (B(p;r),g) of the standard Euclidean 2-spheres § 2 (r) C T p E ~ IR 3 ^. The 

3 The Hessian S = V 2 r is basically the second fundamental form of the immersion U r <— * 
(B(p; r), g). We use the equivalent characterization of shape tensor in order to avoid confusion 
with the standard second fundamental form of use in relativity. 
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shape tensor S measures how the bidimensional metric g^ induced on U r by 
the embedding in (B(p;r),g) rescales as the radial distance r varies. Further 
properties of the equations (|l^) and ( H ) that we need are best seen by us- 
ing polar geodesic coordinates and harmonic coordinates. Recall that normal 
exponential coordinates at p are geometrically defined by 

exp; 1 : B(p; r) -» T p E ~ R 3 (16) 
q i — > exp" 1 (g) = (y 4 ), 

where (y l ) are the Cartesian components of the velocity vector v € T p T, charac- 
terizing the geodesic segment from p to q. Such coordinates are unique up to the 
chosen identification of T P E with M 3 . Since we are dealing with a radial rescal- 
ing, for our purposes a suitable identification is the one associated with the use 
of polar coordinates in T p E ~ R 3 . We therefore introduce an orthonormal frame 
{ei, e2, 63} in T p E such that e\ = d r and with {e2, e^} an orthonormal frame on 
the unit 2-sphere S 2 (l) C T P E. We can extend such vector fields radially to the 
whole Tp£; we consider also the dual coframe {6 2 , 8 3 } associated with {e2, 63}. 
The introduction of such a polar coordinate system in T p E is independent of 
the metric g and thus is ideally suited for discussing geometrically the r-scaling 
properties of (o) B( P - r )- ^ we pull-back to T p S the metric g of B(p; r) C S, we 
get 

g = dr 2 +g(e A ,e B )0 A e B , A, B = 1,2, (17) 

where the components gAB = g{eA, sb) = g^fr, j?, ip) are functions of the polar 
coordinates (§, (f) in T p E ~ K. 3 , associated with the coframe {6 2 ,8 3 }. Note 
that such a local representation of the metric holds throughout the local chart 
(B(p; r), exp" 1 ) and not just at p. In Cartesian coordinates in T p E ~ R 3 one 
recovers the familiar expression^: 

9ab = 5 ab - l -Raku{ P )v k V l + O (\y\ 3 ) . (18) 

In polar geodesic coordinates (r, <p) the equations (|l3|) and (Q take the 
explicit form 

d r S j — S i».S j R r j rJ (19) 
d r gij = 2S\g k j, 

where R l r j r denotes the radial components of the curvature tensor. By taking 
the trace of both equations we get 

d r S = -S\S\ - Ric(d r , d r ), (20) 
gVd r9ij = 25, 

4 Latin indices run through 1, 2, 3; we adopt the summation convention. 
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where Ric(d r , d r ) = R l rir denotes the d r component of the Ricci tensor, and 
S = S l k 8 k i is the rate of volume expansion of g^ (r, ip), (the first equation 
in ( |2C| ) is nothing but the Jacobi operator coming from the second variation of 
the area for g^ 2 '(r, ip)). If we assume that the curvature R(X 7 Y)Z shall then 
be given by fixing the ($, ^-dependence in the factorized metric (p"7|), we can 
consider ( |l9| ) as a system of decoupled ordinary differential equations describing 
the rescaling of the geometry of B(j>\ r) in terms of the one-parameter flow of 
immersions S 2 (r) i— > (U r ,g( 2 " , (r,-&,ip)). Note in particular that the shape tensor 
matrix (S i j) is characterized, from the first set of equations (|l9|), as a functional 
of the given curvature tensor {R 1 ^,^). In such a setting the equation d r gij = 
2<S i9kj can be interpreted by saying that the metric rescales radially along the 
(curvature-dependent) shape tensor S\. In a sufficiently small neighborhood 
of the pole p, from the expression ( |l8| ) and from d r = ^y l c\, we get that in 
Cartesian coordinates we can write 



2 2 
d r g a b = --rR ab (p) - - 



r 3 



Raibk 

(p) + 0(r 2 ), (21) 



where we have evidentiated the trace-free part [y l y — |# r ] of the Cartesian 

tensor y l y k . Roughly speaking, the term | 2^2 i<S jfc r R a ibk(p) represents 

the shear part of the rescaling of the geometry of the spherical surfaces U r == 
{expp(rw) S S : r = const.} as r varies in the neighborhood of the pole p. The 
term —^rR a b(p) is instead responsible for the radial rescaling in the geometry of 
U r . Thus, one may say that in a sufficiently small neighborhood of the point p 
the metric rescales radially in the direction of its Ricci tensor. This latter remark 
will turn out quite useful in understanding the geometric rationale behind the 
choice of a proper averaging procedure for the geometry of (S, g). 

Even if polar geodesic coordinates suggest themselves as the most natural 
labels for the points of B(p; r), they suffer from the basic drawback that their 
domain of definition cannot be a priori estimated and strongly depends on the 
local geometry of (T,,g). In this connection, a much better control on the 
geometry of the balls B(p; r), and hence on (g) B ^ p . r y can be achieved by labelling 
the points exp" 1 ^) £ Be(0,t) with harmonic coordinates, i.e., a coordinate 
system {X 1 } such that the coordinate functions X 1 are harmonic functions 
with respect to the Laplacian on (£,<?). We can do this by starting from the 
given (Cartesian) normal coordinates {y k }, and look for a diffeomorphism on a 
sufficiently small Euclidean ball Be(0, r) C K 3 , 

^:B E (0,r)^B E {0,r) (22) 
y fc = exp; 1 (g)^$(y fe ) = X I 

such that 



^ k =V^ 9 < V ^^ 9ijd ^ k ) = ° (23) 

® k \dB E (0,r) =Id. 
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The standard theory of elliptic partial differential equations implies that the 
harmonic functions so characterized do form a coordinate system in Be(Q,t). 
The important observation is that such harmonic coordinates can be introduced 
on balls of an a priori size as soon as the manifold (E, g) has bounded sectional 
curvature and its injectivity radius is bounded below. Note that from exp" 1 : 
B(p;r) — > TpT, and a generic (i.e., not necessarily harmonic) diffeomorphism 
$ : Be(0, r) — > -Be(0,7") we can define the map 

i exp — 1 (f> 

F^ocxp; 1 : (B(p;r),g) -4 (B E (0,r),S) ^ (B E (0,r),S), (24) 

where S denotes the Euclidean metric on Be(0, r). To any such map we can 
associate the scalar "energy density" JlOf 

e(x) = \\dF\\ 2 = g^S lm ^-^-, (25) 

where the {x 1 } denote given generic coordinates of the points q of (B(p;r), g) 
and where 

is the differential of the map F, thought of as a section of the bundle 

T*X!® F _1 !TR 3 L, , . (27) 

Note that the pull-back bundle F _1 TR 3 over B(p;r) inherits the flat metric 
Si m (F(x)). Thus T*E ® F~ 1 TIR 3 | B(p . r) is naturally endowed with the metric 

g y on T*E and (5; m on F _1 TK 3 , a metric which characterizes the norm ||dF|j 2 . 
Alternatively, e(x) is the trace with respect to the metric of T*E of the pull- 
back by F of the flat metric tensor (Be(0, r), 6). Such remarks emphasize the 
fact that, even if the value of e(x) seems to depend on the choice of local 
coordinates, it is indeed a well-known geometrical quantity: the harmonic map 
energy density associated with the map F : (B(p;r),g) — > (Be(0, r), S). As to 
its meaning in our setting, note that we can always choose the given generic 
coordinates {x 1 } to be the original geodesic coordinates {y 1 }. By exploiting 
such a remark, from the asymptotic expansion (E^), we easily get that in a 
sufficiently small neighborhood of the point p, the scalar e(x) can be written as 

<y) - (> + - 3 R a k b i(p)y k y l + o(\y\ 2 )) Sr* ww - (28) 



If we further assume that the diffeomorphism $ : Be(0, r) — ► Be(0, r) is the 
identity, (i.e., F reduces to exp" 1 ), then = S r a , and we get 

e(2/)lir=e X p- - 3 + l -R kl (p)y k y l + 0(\y\ 2 ) (29) 
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showing that e(x)| F _ exp -i is, roughly speaking, a measure of the curvature 
associated with (B(p; r), g; exp" 1 ). As we shall see, the quantity e(x) will play 
a basic role in understanding the geometric rationale for a proper averaging 
strategy of the geometry of (E,<j>). 



Guided by such geometrical features of geodesic balls, let us go back to the 
study of the scaling properties of (q)b(vt)' T° ^ n ^ s en dj f° r an y r an d s such 
that r + s < injs, let us consider the one-parameter family of diffcomorphisms 
(geodesic ball dilatations) 

H s : (E,p) -» (E,p) (30) 
9 = exp p [r,j((9 r , e 2 , e 3 )] i — ► ff g (g) = exp p [(r g + s)(<9 r , e 2 , e 3 )] 

defined by flowing each point q £ B(p;r) a distance s along the unique radial 
geodesic segment issued at p £ E and passing through q. Let us remark that, 
for any r such that ro < r < injjst, we can formally write 

fl(p;r) =i/ (r _ ro) B(p;r ). (31) 

Thus, in a neighborhood of -B(p; ro) and for sufficiently small r, we get 

M(B(p;r))= f Q d» g = f H( r _ ro) ( Q dn g ), (32) 

JB(mr) JB(v.rn) 



where H* r _ r ^d(gfi g ) is the Riemannian measure obtained by pulling back gdfi g 
under the action of H( r - ro y By differentiating ( |32| ) with respect to r, we have 



—M(B(pr)) = lim 

dr h^o h 



= lim 

h->0 



/s(p;r ) H {r-r a )+h(S d N) ~ Js(p;r ) H (r-r )(^ g ) 



(33) 



Since Hi r _ ro \+h = Hi r - rQ ) ° -ffft> we can write the above expression as 



lim 

h->0 



B(p;r ) 



— lim 



[HhiQdVg) - (Qdfi g 



B(p;r) 



B(p;r) 



lim 



[H*(gdn g ) - (gd/Xg)] 



(34) 
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from which it follows that 



^-M(B(p;r)) = f L dr (gd^ g ) (35) 

ar JB(p;r) 

= ( ( -7T- Q + gdiv(d r )\ dfig 
JB( P -r) KCr ) 



B(p;r) 



where Lg r and div(d r ) denote the Lie derivative along the vector field d r and its 
divergence, respectively, and where we have exploited the well-known expression 
for the derivative of a volume density, 

§- r Vd=lVd9 ab §- r 9a b - (36) 

With these preliminary remarks along the way, it is straightforward to com- 
pute the rate of variation with r of the local average (g) B ( p . r y since (^) implies 

i {Q) ^ = (h) BM (37) 



+ \ { 89ah h«) BM 4 <*W) frl**, Bip . r) 

where (/) B / p . r \ denotes the volume average of / over the ball B(p; r). Explicitly, 
by exploiting ([l9|), we get 

^ (e) B (y,r) = \ j^/ + ^Btor) - (e) B (r,r) ( S ) B(p-r) ■ ( 38 ) 

Thus, the local average (o) B r p . r \ feels the fluctuations in the geometry as we 
vary the scale, fluctuations represented by the shape tensor terms 

(QS)B(p;r) ~ (^)fl B (0;r) B(p;r) (^9) 



governed by the curvature in B(p\r) according to (|20j), and expressing a geo- 
metric non-commutativity between averaging over the ball B(p; r) and rescaling 
its size. Since the curvature varies both in the given B(p\ r) and when we con- 
sider distinct base points p, the above remarks indicate that the local averages 
(f?) B ( p . r ) are subjected to the accidents of the fluctuating geometry of (£,<?). 
In other words, there is no way of obtaining a proper smoothing of g without 
smoothing out at the same time the geometry of (£,<?). 

2.2 Eulerian averaging and Lagrangian scaling 

The use of the exponential mapping in discussing the geometry behind the local 
averages (o) B r p . r ^ makes it clear that we are trying to measure how different 
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the (q) B(p;r) arc lrom tne standard average over Euclidean balls. In so doing 
we think of exp p : T p S — > £ as maps from the fixed space Be(0, r) into the 
manifold (£,<?). In this way we are implicitly trying to transfer information 
from the manifold (£, g) into domains of R 3 which we would like to be, as far as 
possible, to be independent of the accidental geometry of (£, g) itself. Indeed, 
any averaging would be quite difficult to implement if the reference model varies 
with the geometry to be averaged. As already emphasized, this latter task is only 
partially accomplished by the exponential mapping, since the domain over which 
exp p : TpYi — > £ is a diffeomorphism depends on p and on the actual geometry 
of (£, g). We can go a step further in this direction by considering not just the 
given (£, g) but rather a whole family of Riemannian manifolds. To start with 
let us remark that, if we fix the radius tq of the Euclidean ball Be(0; ro) C T p £ 
and consider the family of exponential mappings, exp( p m : T p £ — > (£,<7 (/?)), 
associated with a corresponding one-parameter family of Riemannian metrics 
g a b{(3), < (3 < +00, with g a b{(3 = 0) = g a t, then B(p; ro) becomes a functional 
of the set of Riemannian structures associated with g{(3) a b, < (3 < +00, i.e., 

B(p; r ) 1 — > B p {p- r ) = exp (p/3) [B E (0; r )] . (40) 

In this way, instead of considering just a given geodesic ball B{p; ro), we can 
consider, as (3 varies, a family of geodesic balls Bp(p; r ), all with the same ra- 
dius ro but with distinct inner geometries g a b{(3)- Since J5^ = o(p; ro) = B{p; ro), 
such balls can be thought of as obtained from the given one B(p; ro) by a smooth 
continuous deformation of its original geometry. Under such deformation also 
(g) B(p-ra) becomes a functional, (ff)s (3 ( p ro )j of the one-parameter family of Rie- 
mannian structures associated with g a b{l3), < f3 < +00. The elementary but 
basic observation in order to take properly care of the geometrical fluctuations 
in (q)b(p-to) 1S ^ at tne right member of ( ^7|) has precisely the formal struc- 
ture of the linearization (i.e., of the variation) of the functional {q) B(p- ro ) m ^hc 
direction of the (infinitesimally) deformed Riemannian metric d/d{3[g a b(P)], 
viz., 



df3 ® B t>toro) ~ ( df3 e 



Bf3 (p\r ) 



e 9 ab (P)-^9ab((3)) -(9) B0{p :r o) (9 ab W-^9abW) , (41) 

where the ball i?^(0;ro) is kept fixed while its image B(p;ro) is deformed ac- 
cording to the flow of metrics g a b(0), < (3 < oo. 

In a rather obvious sense, (^) represents the active interpretation corre- 
sponding to the Eulerian passive view associated with the ball variation B(p; r ) — 
B(p;r). In other words, we are here dealing with the Lagrangian point of 
view of following a fluid domain in its deformation, where the fluid particles 
here are the points of B(p; ro) suitably labelled. This latter remark suggests 
that in order to optimize the local averaging procedure associated with the lo- 
cal average (g)B( P -r)' instead of studying its scaling behavior as r increases, 
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and consequently be subjected to the accidents of the fluctuating geometry of 
(B(p; r), g a b), we may keep fixed the domain B(p; ro) and rescale the geometry 
inside B(p; r ) according to a suitable flow of metrics g a b — * g a b{(3), < (3 < oo. 
Correspondingly, also the matter density g will be forced to rescale g — > g((3), 
and if we are able to choose the flow g a b — > g a b{/3) in such a way that the lo- 
cal inhomogeneities of the original geometry of (E, g) are smoothly eliminated, 
then the local average (g) B ( p . r ) comes closer and closer to represent a matter 
averaging over a homogeneous geometry. 

In order to put such Lagrangian picture on a firmer ground we need to 
consider the important issue of how to label the points of (B(p; r ), g{(3)) during 
the deformation process. For (3 = the natural choice are (polar) geodesic 
coordinates associated with exp p . However, as the metric changes along the 
flow (B(p;r Q ),g(f3)), the remarks of the previous section suggest that polar 
geodesic coordinates may not be the optimal choice as (3 varies, since we do not 
have an a priori control on the domain of validity of such a coordinate system. 
However, if we select the flow g a b — > g a b{(3), < (3 < oo, in such a way that the 
sectional curvatures of g{(3) stay bounded and the injectivity radius is bounded 
below as (3 varies, then such an a priori control can be achieved by labelling 
the points of expj" 1 ^ by harmonic coordinates. Thus, in order to properly label 
the points of (B(p;r ), g{(3)), instead of considering the one-parameter family 
of maps expj" 1 ^ : (B(p; ro), g((3)) — > (-Be(0; ro), S), we shall consider the family 
of maps 

F(f3) = $(/?) o exp^ : (B(p; r ), <?(/?)) ^ (B E (0; ro), 5) ^ (B E (0; r ), S), 

(42) 

where $(/3) : (Be(0; r ), 5) — > (Be(0; ro),S) is a /3-dependent diffeomorphism on 
(Be(Q; ro), 5), (inducing a coordinate change in (-Be(0; ro), 5)). By mimicking 
(p5|), we associate the harmonic map energy density with the family of maps 
F{(3) providing the Lagrangian coordinates that allow to follow the points of 
(B(p;ro),g(/3)) during the deformation, i.e., 

e(J3;x) = \\dF((3)\\ 2 = {P^ J*^ } , (43) 

where the {x 1 } denote the fixed generic (Eulerian) coordinates of the points of 
(E p , g{(3)). If we compute the rate of variation with (3 of e(/3; x) we get 

Since we eventually want to label the points of (B(p;ro), g{(3)) with harmonic 
coordinates, the most natural thing to do is to deform F((3) according to the 
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harmonic map flow, i.e., 



where 



^F m (P) = ApF m ((3) 
F m {P = 0) = exp" 1 , 



(45) 



AgF m (/3) = g ij (P) 



(46) 



is the harmonic map Laplacian associated with the Riemannian connection 
F k ij(g(P)) of (B(p;r ),g(p)). In this way F m {j3) evolves towards harmonic 
maps (i.e., the Lagrangian coordinates labelling the points of (B(p;r ),g(P)) 
will be harmonic coordinates as (3 — > oo). Note that ([|5]) is a strictly parabolic 
initial value problem admitting a solution for sufficiently small j3 > 0, (in our 
case, standard theorems in harmonic map theory imply that the solution actu- 
ally exists for all < (3 < oo). If wc insert (45) into (B) we get 



d_ 

dp~ 



e((3;x) 



dF l {(3)dF m {(3) d 



+2g ij (/3)5 l r 



da* 
dF l (J3) d 
dx % dxi 



dxi dp 



(47) 



(A F m ((3)). 



Now, note that a Bochner type formula allows to compute the Laplacian of 
e(/3; x) according to 



A 



(e(p;x)) = A (\VF(p)f 



= 2 \W 2 F(/3)\ Z + 2Ric(VF(p),VF(p)) + 2g (VF(p), VApFtf)) 
where the notation is as follows: 

|V 2 ^(/3)| 2 = g ik (p)g^ l (P)S rs V 2 i F r (P)Vl l F s (P), 



%F r {(3) 



Ox 



g(VF(P),VA F(P)) ± g* k (P)5 rs V(A F r (P))\7 k F s {P), 
Ric(VF(p),VF(p)) = 6 lm dF y ] ] R^{(3). 



dx % 



dxi 



Equation mM implies that we can write (]47[) as 



d_ 

dp 



e(P;x) 



dF l (f3)dF m {(3) 8 



-2 \V 2 F{(3) 



dx l 

,2 



26. 



dxi dp 
dF l {P) dF m (P) 



lm ~ 



g v(P)+A (a(/3;x)) 
R l i{P). 



(48) 

(49) 
(50) 
(51) 
(52) 

(53) 



dxi 



This latter expression shows that, as expected, the /3-evolution of e(/3; x) is 
strongly influenced by the chosen evolution of g % i(P). However, such an evo- 
lution becomes particularly simple if we choose the flow of metrics g lJ (P) in a 
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way that is strongly reminiscent of the radial scaling properties of the metric 
©, viz., 

^g i i([3)=2RV(l3), (54) 
or in terms of the covariant components gij {(3) 

^9iM = -2i2yC3). (55) 
Under such a choice (|5^) reduces to 

^e(/3; x) = A p (e(/3; a,)) - 2 |V 2 i^(/3)| 2 , (56) 

which implies that, as (3 increases, the maximum of e(/3;x) is (weakly) decreas- 
ing. According to the geometrical meaning of e{(3;x), (see (p9|)), curvature 
inhomogeneities tend to decrease along the flow (^35|), at least in a sufficiently 
small neighborhood (in the C fc -topology) of (B(p;ro), g(/3)). 

Obviously, these remarks are strongly reminiscent of the properties of the 
Ricci flow on a Riemannian manifold (£,<?), i.e., 

£j9ab{ri) = -2Rab{ri) ^ 

9ab(V = 0) = 9ab, 

studied by Richard Hamilton and his co-workers in connection with an analytic 
attempt to proving Thurston's geometrization conjecture. As is well-known, 
the flow (|57| ) is weakly-parabolic, and it is always solvable for sufficiently small 
■q. Moreover, it is such that any symmetries of g a b{v = 0) = g a b are preserved 
along the flow. The flow (|57|) may be reparametrized r\ — > A by an ?y-dependent 
rescaling and by an ^-dependent homothety g a b{g) — * ffafc(A) so as to preserve 
the original volume of (£,<?). In this latter case in place of ( p7| ) we get the 
associated flow 

f M9abW = -2R ab (X) + § (£(A)) e(a) 9abW 
I 3ah(A = 0) = g ab , 

whose global solutions (if attained) 

9ab = , lim 5afc(A) (59) 

A — >+oo 

are constant curvature metrics. 

Thus, a Lagrangian picture shows in a very direct way that the natural flow 
gab — >• g a b(P) providing the simultaneous rescaling of the matter distribution and 
the geometry is a (suitably normalized) Ricci flow. Clearly, its normalization 
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cannot be selected simply on the basis of geometric convenience {e.g., fixed 
volume, as in the standard Ricci-Hamilton flow(|5^)), and we can realize our 
matter averaging program only if there exists a global solution g a b — > g a b(/3), < 
f3 < oo of the Ricci flow which preserves the total matter content M(S) of (E, g). 
Such a problem together with the set up of a renormalized effective dynamics of 
smoothed out cosmological spacetimes is the headline of our forthcoming work. 
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